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ABSTRACT 

The high quality spectra required for radial velocity planet searches are well- 
suited to providing abundances for a wide array of elements in large samples 
of stars. Abundance ratios of the most common elements relative to Fe are 
observed to vary by more than a factor of two in planet host candidates. This 
level of variation has a substantial impact on the evolution of the host star and 
the extent of its habitable zone. We present stellar models of 1M Q stars with 
custom compositions representing the full range of these non-solar abundance 
ratios. We find that the effects derived from variation over the observed range of 
[O/Fe] has a particularly dramatic effect. Habitability lifetimes for some classes 
of orbits can vary by gigayears for the observed range in [O/Fe]. 

Subject headings: stars: abundances — stars: evolution 
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Introduction 



In general stellar models are created for different total metallicities scaled relative to 
the solar abundance pattern, but variations in the abundance ratios at a given [Fe/H] are 
not considered, except for uniform enhancement of the a elements (O, Ne, Mg, Si, S, Ar, 
Ca, Ti) at very low metallicities. This is an important problem in its own right, but recent 
developments in the field of extrasolar planets have added another layer of interest. 



Recent discoveries by Kepler o 



Charpinet et al. 



2011 



Fressin et al 



zones of stars (i.e . Kepler-22 b (IBorucki et al 



Eart h-sized and smaller planets (iBorucki et al. 



2011 



2012) and planet s near or within th e nominal habita 



Die 



Mayor et al. 



20091 ) GJ 667C c and d flDelfosse et al 



2012 



, GL 581 c and d flUdrv et al 



2012 



Anglada-Escude 



2007 



20121)) lends 



a sense of immediacy to the quest ion of what makes a planet habitable . The present 



position of a star's habitable zone ( jSelsis et al. 



2007 



Lammer et al. 



20091 ) can be estimated 



dire ctly from observations and appropriate assumptions about planetary atmosphere s 



see 



Wordsworth et al 



2011 



Hu fc Ding 



2011 



Kaltenegger. Segura. Sz Mohanty 



2011 



for 



atmospheric modeling of GL 581 d). The evolution of habitable zones over time is also 
of fundamental importance. Although life arose quickly on Earth, other steps of critical 
astrobiological interest took considerably longer. Potentially detectable atmospheric changes 
produced by the biota, in particular the presence of free oxygen, methane, or nitrous oxide, 
occurred after roughly a billion years. Complex life required a further three billion years. 
Arguments have been proposed that suggest certain evolutionary breakthroughs are low 
proba bility events, making it likely that evolu tionary timescales will be long as a general 



rule. (ILineweaver fc Davis 



2002 



Carter 



2008 



and [Fe/H] is a discipline of long standing (IMowlavi et al. 



). The variation of stellar lifetimes with mass 



2012 



for a recent example) 



and changes in 



Muirhead et al. 



labitable zone extents with mass and [Fe/H] have been considered (e.g. 



20121 ). but other abundances may have a substantial effect. 
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High quality spectra required for radial velocity planet searches are well-suited for 
determining stellar abundances. For large samples of stars it is possible to statistically 
separate the i ntrinsic variation of elemental abundanc e ratios from observational 



uncertainties ( iDe Silva et al. 



2006 



Stuart fe Ord 



20091 ) . An analysis of published 



abundances 



( jBond et al. 



or thirteen elements derived from the AAT radial velocity planet search 



2006 



20081 ) provides measurements of the range of variation of abundance 
ratios relative to Fe for a hundred stars. The sample consists of mostly G and some 
F and K dwarfs with ages > 3 Gyr and -0.5<[Fe/H]<0.3, providing an observationally 
uniform sample that approximates as well as possible the primordial stellar abundances. 
The elemental ratios relative to the sun ([X/Fe]) vary substantially, more than can be 
accounted for by observational error. The 3ui n t r insic variation (this spread is an rms error, 
not a standard deviation for a normal distribution) is close to a factor of two about the 
mean for C, O, Na, Mg, Al, Si, Ca, and Ti, from 1.7 for Na to 2.3 for O. Si, Cr, and Ni 
have smaller spreads of a factor of 1.4-1.5 about the mean. The r- and s-process elements 
(species significantly heavier than Fe produced in evolved stars and supernovae) have even 
larger variations. Also, while these ranges encompass the majority of nearb y planet host 



2008|). The star 



candidates, there are even more extreme outlier systems (jBond et al. 
HD53705 has [O/Fe] = 4, more than ■;>//,„/,.,„.■,„• higher than the mean. HD136352 is similarly 



enhanced in [C/Fe] (jBond et al 



20061 ) . Both of these are low metallicity systems ([Fe/H = 
-0.18 and -0.26, respectively). Enhancements of this magnitude will substantially increase 
the range of orbits that remain habitable for a multi-Gyr window. 

Stellar evolution calculations are sensitive to the assumed composition. This sensitivity 
arises from two c omposition-dependent effects: the equation of state (EOS ) and the 



radiation opacity ( [Rogers. Swenson. fe Iglesias 



1996 



Iglesias fe Rogers 



19961 ). The equation 



of state is primarily affected by the number of free particles (less when more nucleons 
are in heavy nuclei) and the effects of electron screening, where the mutual repulsion by 
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the positive charge of high Z nuclei is reduced by intervening electrons. The changes 
in the EOS are minor, and much less important than the opacity. The elements C, N, 
O, Ne, Si, Ti, and Fe-peak (stable elements with atomic numbers close to Fe) elements 
are particularly important because of their relative abundance and number of electron 
transitions. Rearranging the proportions of different species at a co nstant metallicity can 



result in opacity changes of tens of percent (jlglesias fe Rogers 



19961 ) . Opacity changes the 



rate of leakage of radiation. Increased radiation pressure in the stellar envelope drives 
expansion, resulting in larger radii and lower effective temperatures. Also, a slower rate of 
energy loss requires a slower rate of nuclear burning to maintain hydrostatic equilibrium. So 
we expect stars with enhanced abundances to be cooler, lower luminosity, and longer lived. 
The large spread in O is particularly interesting. Given that its high abundance causes 
it to play a role second only to Fe in stellar opacity, stellar models with large variations 
in O abundances might be expected to show a significant impact on stellar evolution and 
observable quantities. 

We explore whether the observed level of abundance ratio variations cause significant 
changes in the evolution of stars and therefore their habitable zones. Section [2] describes 
the simulations used, Section [3] presents the quantitative changes in the stellar evolution 
and habitable zones, and Section H] examines implications of the work. 



2. Simulations 



For this w ork one solar mass sta rs were simulated using the stellar evolution 



code TYCHO ( lYoung fc Arnett 



20051 ). TYCHO is a ID stellar evolution code with a 



hydrodynamic formulation of the stellar evolution equations. It uses OPAL op acities 



(jlglesias fe Rogers 



1996 



Alexander fe Ferguson 



1994 



Rogers fe Navfonov 



combined OPAL and Timmes equation of state (HELMHOLTZ) (ITimmes fe Arnett 



2002h. a 



1999; 
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Rogers fc Nayfonovll2002l ). gravitational settling (diffusion) (IThoul. Bahcall. fc Loeblll994f ). 
general relativistic gravity, automatic rezoning, and an adaptable nuclear reaction network 
with a sparse solver. A 177 element network terminating at 74 Ge is used throughout 



the evolution. The network uses the 



2000 



Angulo et al. 



1999 



Iliadis et 



Langanke fc Martmez-Pinedo 



al 



2001 



Wiescher et al 



atest REACLIB rates (Rauscher fc Thielemann 



200 6]), weak rates from 



Graboske et al. 



( jl973h . Neutrino 



2000), and screening from 
cooling from plasma processes and the Urea process is included. Mass loss is included but is 
trivial for a lM ff m ain sequence star. (Heightened early mass loss seen in some young stars 



( Wood et al. 



( IMeakin &: Arnett 



2005D is not included.) It incorporates a descript i on of turbulent convect ion 



2007 



Arnett. Meakin. &: Young 



2009 



2010 



Arnett fc Meakin 



20111 ). It 



has no free convective parameters to adjust, unlike mixing-length theory. 



A solar composition from 



Lodders et al. 



( 120091 ) was adjusted to the mean abundance 



ratios of the AAT sample for the elements observed, with all other elements being 
maintained at solar values for the standard composition. The adjusted values are listed 



in Table [TJ The i ntrinsic physica 



are accounted for ( iDe Silva et al. 



spreads u 



2006 



intrinsic. 



re maini ng when observational errors 



Stuart fc Ord 



20091 ) in the element to iron ratios 



are the basis for the varian t models. The methodology is described fully in a forthcoming 



paper ( jPagano et al. 



20121 ) . but the relevant values are quoted here. Table [Tj gives the 
enhancement /depletion in X/Fe for each element. In each model one element was varied to 
the 3iii n trinsic extremes of the distribution, giving one standard composition model and two 
models for each of the elements C, O, Na, Mg, Al, Si, Ca, Ti. Of the set of elements in the 
observations mentioned in Section [TJ we do not model variations in Cr and Ni because of 
their very small ranges or the r- and s-process because of their very low abundance. We 
do include Si in spite of its small variation because it has a relatively high abundance. All 



other elements are held constant. 
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Determinations of O abundan ce are notoriously sensitive to non-LTE effects 



( IGrevesse. Asplund. fc Sauval 



20071 ). which may introduce random unphysical errors in the 
O abundance measurement not accounted for in the quoted observational errors, so we 
also examine a smaller range of O variation. Inste ad of choosing a random value, we use 



Ramirez. Melendez. fc Asplundl (120091 ) . 



the Ui n trinsic for a sample of 40 "solar twins" from 
Qualifying as solar twins requires observable parameters quantitatively close to the sun. 
In order for a star to match the sun closely in physical observables, it is necessary for it 
to have a composition close to solar. Using Ui ntr insic for the solar twin sample therefore 
provides a very conservative estimate of [O/Fe] variation. 

Custom opacity ta 



( llglesias &: Rogers 



:>les for each composition were obtained from the OPAL project 
19961 ). All models were run until core hydrogen exhaustion. The models 
provide high time resolution histories of luminosity and effective temperature. Translation 
of these quantities into the extent of a habitable zone is anything but straightforward. 
Habitability depends on, among other factors, atmospheric composition, albedo, density 
structure and therefore planet mass, cloud cover, assumptions about water loss, geophysical 
exchange with the atmosphere, and assumptions about runaway greenhouse effects. For 
this work we chose a reasonably conservative prescription for an Earth mass planet that 
avoids the most extreme predictions for boundaries. We chose to estimate the inner and 



outer r adii of the habitab le zone using the prescription of 



(11993h 



Selsis et al. 



Kasting. Whitmire. fc Reynolds 



(120071 ) . The radii depend upon stellar T e ff and luminosity, with 



and 



lout (JoutQ ^outX* ^out^-n. 



L \ 5 



(2) 



where a in = 2.7619 xlO" 5 , b in = 3.8095 xlO" 9 , a out = 1.3786 xl0~ 4 , b out = 1.4286 xlO" 9 , 
and T* = T e ff — 5700. For a 0% cloud cover atmosphere and assuming an inner boundary set 
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by the water loss limit, l in& = 0.95AU. For a 50% cloud cover atmosphere, l ou t& = 1.95AU. 
We choose the larger inner boundary in order to give a 1 AU orbit for a solar composition 
a habitable lifetime of 5.5-6 Gy. With this type of treatment, the habitable zone edges 
can be easily scaled from the stellar quantities for different assumptions, and the relative 
differences between the stellar models are robust. 



3. Results 

The impact of abundance variations for several elements are large enough to be 
easily seen in the HR diagram, in which the evolutionary tracks of the stellar models are 
plotted in log L/L & versus effective temperature T e ff(K). Tracks are shown in Figured] 
for enhancements (left) and depletions (right) in C, O, Na, Mg. In all cases enhanced 
abundance ratios cause larger changes than depletions. C, Na, and Mg have small but 
noticeable effects. C has a high abundance but relatively few electron transitions and low 
ionization potential. Mg and Na have lower abundances but higher opacity per gram than 
C, resulting in a similar degree of shift in the tracks. Si has less impact due to its smaller 
range of variation, and Al, Ca, and Ti have very small effects due to their small abundances 
and can be neglected for our purposes. Although the main sequence tracks are similar 
in morphology, substantial differences arise on the pre-main sequence (pre-MS). This is 
because the changed C and O abundances affect the energy release and therefore transient 
convective core size, during the partial CN cycle burning that occurs during the pre-MS 
phase. This stage is, however, too short-lived to be relevant to habitability. 

The largest changes, unsurprisingly, arise from variation in O. Figure [2] shows the 
evolutionary tracks for the Standard, enhanced, moderately enhanced and depleted, and 
depleted O compositions. As with the other elements, enhancements result in larger changes 
than depletions due to the larger absolute changes in mass fraction of the element. For 
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this letter we will concentrate our analysis on the adjusted O abundance ratios since the 
evolutionary tracks for O represent the extremes of variation. 

Both the luminosity and T e ff of enriched compositions are systematically lower at a 
given age. The peak temperatures reached on the main sequence are \ogT e ff = 3.753 for 
enhanced [O/Fe] and logT e jj = 3.78 for low [O/Fe] (logT e jj = 3.762 for the sun). This 
gives a temperature difference of ~360K, or a 6% deviation from solar effective temperature. 
The temperature range for the moderate O cases is 135K, or 2.3%. 

The time evolutions show the most profound effect to be on the pace of the evolution. 
The main sequence turnoff for the low [O/Fe] model occurs at an age of 9Gy. The standard 
composition (which is somewhat more O-rich than solar) has a turn-off age of 10.5Gy, and 
the O-rich model turns off at 12Gy at lower logL and T e fj- A 33% increase in stellar 
lifetime is interesting in itself, but the effects on habitability require further analysis. 
Orbits with semi-major axes near the edges of habitable zones will not remain there for 
the entire lifetime of the star. Earth itself is in this predicament, with many models 
predicting only about a billion years more of habitability. Figure El shows the inner and 
outer edges of habit able zones for the h igh, low, and standard [O/Fe] as a function of time, 



calculated using the 



Selsis et al. 



(120071 ) prescription. The gray bar marks a 1AU orbit. We 
consider several properties of habitable orbits. First, we find the innermost and outermost 
continuously habitable orbits and innermost and outermost orbits habitable for 4Gy. An 
Earth-centric measure of habitable zone persistence is the habitable lifetime of a 1AU orbit. 
The radial extent of the habitable zone varies substantially over the life of the star as well, 
increasing monotonically. The evolution of the outer edge of the habitable zone in particular 
is steepest and largest in absolute terms for the O depleted model. These parameters are 
provided in Table [2] 
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discussion 



Potential for habitability will be an important co nsideration when allocating resources 



for t 



r e very difficult direct observations of exoplane ts (ITessenyi et al 



2011 



Rauer et al. 



2011 



Kaltenegger fc Traub 



2012 



von Paris et al 



20091 ). A first cut of candidates can be made 



on the current habitability of a planet estimated directly from the measured properties of 
the host star. The current habitability of a planet is not the only consideration in finding 
Earth-like worlds, however. If the ultimate goal is the discovery of life, the evolution of 
habitable zones must be taken into account. Life did not modify the Earth's atmospheric 
composition at levels that might be detectable by near-future astronomical observations 
until about a billion years after its formation. Complex life on Earth developed during the 
late Ediacaran period, 635 - 542 Mya, some 4Gy after its formation. For the holy grail of 
Earth-like complex life, if Earth is at all representative, we should look for planets that 
have spent several billion years inside their habitable zones. 

For a star otherwise identical to the sun in mass and Fe abundance at the lower end of 
the range of [O/Fe] observed in planet host candidates, an Earth-like planet at 1AU from 
the star would not (given the assumptions about planetary atmospheres herein) remain 
habitable long enough for complex life to have developed at the pace that occurred on 
Earth. The effect is obviously dramatic for the extremes of O variation, but there is even 
an effect that could be important to the less extreme systems with planets on the edges 
of their nominal habitable zones. The 1 AU habitability lifetime differs by nearly a billion 
years between the moderately enhanced and depleted cases. This is the variability deduced 
for a set of stars deliberately chosen to be maximally sun-like. The magnitudes of the effect 
for the observed ranges of [C/Fe], [Na/Fe], and [Mg/Fe] are hundreds of millions to a billion 
years. 



Any orbit that is not continuously habitable varies in its habitable lifetime by billions 
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of years between stars of different abundance ratios. If we limit ourselves to systems that 
are continuously habitable for the entire main sequence lifetime of the star for all [0 /Fe] we 
limit ourselves exceedingly. For this atmospheric prescription the only orbits habitable for 
all of the models with the observed range of [O/Fe] at solar [Fe/H] extend from 1.35AU 
to 1.55AU. Ranges are similarly narrow though shifted in radius for other [Fe/H]. For any 
given single composition the range is much larger, increasing for lower abundance ratios. We 
also find the orbits that should have at least 4Gy of continuous habitability. This broadens 
the range of allowed orbits considerably. High [O/Fe] models have the widest range of 4Gy 
habitable orbits by a small margin due to their longer overall main sequence lifetime and 
shallower rate of change of the zone boundaries. A point that is not commonly considered 
is that the majority of these orbits are habitable during the latter part of the star's lifetime, 
meaning that many of the planets that may be favorable to the development of life that we 
discover will have only begun their habitable lifetimes. Also, while these ranges encompass 
the majority of nearby planet host candidates, there are even more extreme outlier systems. 
The stars HD53705 ([O/Fe] = 4) and HD 136352 ([C/Fe] also ~4)have a substantially larger 
range of orbits that remain habitable for a multi-Gyr window with such large enhancements. 

Conversely, this work also reiterates the much greater difficulty in determining the 
extent of a star's habitable zone due to planetary conditions. Our discussion is limited to 
the effects on the star. The pace of a star's evolution can cause a very large change in the 
habitable lifetime of an orbit even if the radial distance of the habitable zone from the 
star changes by only a small amount. At 6.1Gy the habitability of the 1AU orbit ends 
for the standard composition. The inner edge of the habitable zone for the moderately 
enhanced O case at 6.1Gy is at 0.97AU. This 3% difference results in an additional 0.6Gy 
of habitability due to the slower increase in luminosity of the more O-rich composition. 
A very minor change in the model parameters chosen even from the single atmospheric 
prescription we examine would change the habitable zone radius more than this. This 
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also assumes that prescriptions for planetary habitability based on assumptions about the 



evolution o 



Bond et al. 



plane tary atmospheres are appropriate for different planetary compositions. 



( 120081 ) show that the mineralogy of solid bodies arising from protoplanetary 
disks of non-solar abundance ratios can be dramatically different from the solar system. The 
composition of planetary atmospheres and the geological cycles that modify and maintain 
them can be expected to vary similarly. 

The viability of exoplanets at the edges of their habitable zones depends directly on 
the abundances of all the common elements. The stellar characteristics that determine the 
size and persistence of habitable zones, completely apart from properties of the planets 
themselves, change at an important level. Proper characterization of a system requires 
more data than just [Fe/H]. A concerted effort to obtain high quality spectral observations 
of all planet hosts of astrobiological interest would be valuable. 

Acknowledgements: We would like to thank Eric Mamajek and Dave Arnett for 
useful discussions, Mark Richardson, Matt Mechtley, and Vahid Golkhu for abundance 
tables, and the referee for most excellent suggestions for improving the paper. This work 
was supported in part by a NASA Astrobiology Institute grant. 
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Table 1. Abundance Enhancements/Depletions Relative to Standard Composition 



Element 


Depletion 


Enhancement 


LX./H Qfnnrlnrri^ 


C 


0.48 


2.09 


0.173 





0.44 


2.28 


0.051 


Moderate O b 


0.67 


1.48 




Na 


0.58 


1.74 


0.068 


Mg 


0.43 


2.32 


-0.029 


Al 


0.55 


1.81 


-0.124 


Si 


0.71 


1.41 


0.078 


Ca 


0.54 


1.86 


0.027 


Ti 


0.40 


2.49 


0.072 



Standard composition values relative to solar 



b Uses the smaller range of O variation observed in twins 
(see Section [2]). 
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Fig. 1. — Evolutionary tracks for depleted (left) and enhanced (right) abundance ratios 



Black is the standard composition of 



Pagano et al 



(120 121 ) . Na (blue), Mg (cyan), C (red), and 



O (green) are each varied according to Table [T] while the rest of the metals are held constant. 
The effects of abundance enhancements are consistently larger than those of depletions, 
driving stars to lower T e ff and luminosity. O causes the most pronounced change. Kinks in 
the tracks are slight variations in the models' envelope solutions and represent a temperature 
variation of order 5K. 

Table 2. Habitable Zone Extent 



Composition Continuously 4Gy Habitable Width at Width at H Habitable Lifetime 
Habitable ZAMS Exhaustion of 1AU Orbit 

AU AU AU AU Gy 



High 


1.15 - 


1.55 


0.84 - 


■ 2.11 


1.3 


1.7 


9.0 


Medium High 


1.22 - 


1.64 


0.9 - 


2.04 


1.34 


1.79 


6.7 


Standard 


1.25 - 


■ 1.65 


0.92 - 


■ 2.06 


1.35 


1.8 


6.2 


Medium Low 


1.25 - 


1.67 


0.94 - 


■ 2.19 


1.35 


1.8 


5.8 


Low 


1.35 - 


■ 2.12 


1.03 


- 2.2 


1.41 


1.89 


3.5 
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Fig. 2. — Evolutionary tracks for Standard (black), enhanced O (cyan), moderately en- 
hanced O (blue), moderately depleted O (green), and depleted O (red). Similarly to the 
tracks for other elements, luminosity and temperature shifts are smaller for depletions than 
enhancements. 
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0.5 1 1 — ' 1 >— " — ' ■— 

5xl0 9 10 

t(Gy) 



20071 ) versus time for 



Fig. 3. — Inner and outer edges of the habitable zone (jSelsis et al. 
O. Habitable zones are shown for models with depleted (dotted), moderately depleted (short 
dash), standard (solid), moderately enhanced (long dash) and enhanced (dot-dash) O abun- 
dance ratios. Depletions systematically move the habitable zones outward due primarily to 
increased luminosity and secondarily to increased T e ff. The difference in location increases 
markedly with age. The most dramatic effect results from the substantial difference in main 
sequence lifetimes. 
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